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We investigate experimentally and theoretically the third harmonic generated by a double-layer 
fishnet metamaterial. To unambiguously disclose most notably the influence of the magnetic res- 
onance, the generated third harmonic was measured as a function of the angle of incidence. It is 
shown experimentally and numerically that when the magnetic resonance is excited by pump beam, 
the angular dependence of the third harmonic signal has a local maximum at an incidence angle of 
~ 20°. This maximum is shown to be a fingerprint of the antisymmetric distribution of currents 
in the gold layers. An analytical model based on the nonlinear dynamics of the electrons inside 
the gold shows excellent agreement with experimental and numerical results. This clearly indicates 
the difference in the third harmonic angular pattern at electric and magnetic resonances of the 
metamaterial. 

PACS numbers: 81.05.Xj, 42.65.Ky, 73.20.Mf, 42.70.Mp 



I. INTRODUCTION 

The emergent field of metamaterials has brought opti- 
cal materials to a qualitatively new level. It became pos- 
sible to access new functionalities and optical properties 
of media by applying subwavelength structuring. Spec- 
tral selectivity and extraordinary optical transmission, 1 
chirality, 2 anisotropy^ optical magnetism and negative 
refraction could be assigned to thin-film media by means 
of modern nanolithography. 5 9 The majority of effects 
in optical metamaterials arise due to the excitation of 
plasmon polaritons at interfaces between metal inclusions 
and surrounding dielectrics. Plasmon polaritons produce 
highly localized electromagnetic field densities, making 
metamaterials attractive from the point of view of nonlin- 
ear optical effects enhancement 10 13 and tailoring of non- 
linear optical properties. 14 22 Magnetic metamaterials, 
i.e. metamaterials that mimic optical magnetism by sup- 
porting circular current plasmonic modes are of spe- 
cial interest since the respective circular currents can play 
a significant role in the nonlinear optical response ] n i 24 ~ — 
Although considerable effort was directed towards de- 
termining peculiarities of the nonlinear-optical response 
caused by the excitation of magnetic resonances j 10 i 1 1 » 24 

In this paper we attempt to demonstrate an implicit 
evidence of the symmetry-induced characteristics of the 
nonlinear response in magnetic metamaterials. There- 
fore we clearly disclose the magnetic mode contribution 
to the third-order nonlinearity of the fishnet metamate- 
rial. This is done by means of angular spectroscopy of 
third harmonic generation (THG) and numerical model- 
ing of THG with a Fourier modal method (FMM) with 
a nonlinear extension. The results are specific to the 
case when the magnetic resonance of the metamaterial is 



excited by pump radiation. The magnetic mode contri- 
bution arises from the antisymmetric current distribution 
in the two gold layers of the metamaterial and is revealed 
as a local maximum of THG intensity in the angular de- 
pendence at tilted incidence. The data is supported by 
an analytical model based on the dynamics of coupled 
nonlinear oscillators. This reveals the strong influence of 
the resonance symmetry on the third harmonic angular 
radiation pattern because of the retardation effects. 



II. SAMPLE 

The fishnet structure is laterally defined by electron 
beam lithography (Vistec SB350OS) and a lift-off tech- 
nique on a Si02 substrate and comprises a set of rectan- 
gular holes fabricated in a three-layer Au-MgO-Au het- 




FIG. 1. (a) A scanning electron microscope (SEM) pic- 
ture taken from the top view. The inset (b) shows a mag- 
nified SEM picture. The period p is 500 nm in both lat- 
eral directions, (c) The details of geometry parameters are 
w x — 110 nm, w y — 290 nm, d = 23 nm, s = 65 nm. 
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erostructure. The resulting structure has thin wires with 
a width of w x — 110 nm and broad wires with a width of 
w y = 290 nm. The structure has a period of p = 500 nm 
in both lateral directions. The thicknesses of both Au 
films are d — 23 nm and the thickness of the interme- 
diate dielectric MgO film is s = 65 nm. Figure QJa) 
and Fig. DJb) show different scales of a scanning elec- 
tron microscope (SEM) image of the sample taken from 
the top view. The parameters of the sample are shown in 
Fig. Die). The specific design parameters were chosen to 
match the magnetic resonance wavelength to the telecom 
wavelength range. 



III. LINEAR OPTICAL RESPONSE 

Simulated and measured linear absorption spectra of 
the sample for different angles of incidence are shown 
in Fig. [2ja) and Fig. [2jb). The angles vary from 10° 



up to 50° with 10° steps. The linear absorption A\ is 
measured with an integrating sphere module of a Perkin 
Elmer Lambda 950 spectrometer in the spectral range 
of 1.20 /im- 1.80 /im in 2nm steps. For the theoretical 
description of the problem we applied the FMM which 
allows solving the linear diffraction problem for an arbi- 
trary anisotropic bi-periodic multilayer structured For 
the simulation we used the parameters of the sample as 
measured using the SEM. The angular spectroscopy of 
absorption measured with p-polarized incoming light is in 
agreement with the theoretical predictions, i.e. existence 
of an absorption peak at the wavelength of approximately 
A = 1.54 /im for normal incidence which is blue-shifted 
as the angle of incidence is increased. The magnetic mo- 
ment of the resonance results from the currents inside the 
broader wires of the metamaterial flowing in the opposite 
directions. This is demonstrated in Fig. [2fc) where the 
calculated phase difference A</> between the electric field 
in top and bottom gold layers is shown as a greyscale 
plot from the top view. 




y(jum) 

FIG. 2. Numerically simulated (a) and measured (b) linear 
absorption A\ as a function of wavelength A in the spectral 
range of the magnetic resonance, plotted for different angles 
of incidence: 6 = 10°, 20°, 30°, 40° and 50°. The dotted line 
in the simulation plot shows the absorption for normal inci- 
dence, (c) The calculated phase difference between the 
^/-components of the electric fields in top and bottom gold lay- 
ers at the center of the resonance peak for normal incidence, 
at 1.54 /im. 



IV. NONLINEAR OPTICAL RESPONSE 

For the nonlinear measurements a setup based on an 
optical parametric amplifier (OPA) was used operating 
at wavelengths of 1.49 /im, 1.54 /im, 1.56 /im and 1.60 jam 
and having an average output power of 3 mW focused to 
a 300 /im-spot from the air side of the sample. The OPA 




FIG. 3. The setup for angular spectroscopy of the third har- 
monic generation (THG) intensity. The pump polarization is 
set to p and p-polarized third harmonic radiation is detected 
with a photomultiplier tube (PMT). The diffraction in the x 
direction is not shown. 
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FIG. 4. Subplots (a)-(d) show the third harmonic signal as 
a function of the angle of incidence for different wavelengths 
in the spectral vicinity to the magnetic resonance. For com- 
parison (e)-(h) show the linear absorption A\ at the same 
fundamental wavelengths and (i)-(l) shows the linear trans- 
mission T at the corresponding third harmonic wavelengths. 
The vertical dashed lines indicate the angular positions of 
the appearance and the disappearance of diffraction orders. 
The black dots represent the experimental data and the dot- 
ted lines represent the simulation results. The solid lines are 
curves calculated with Eq. (fl"2]) . This equation represents an 
analytical model which describes the nonlinear response of 
coupled oscillators. 



was pumped by a Nd:YAG laser with pulse duration of 
5ps and a repetition rate of 5 kHz. The resulting fluence 
took values up to 700 /iJ/cm 2 in the plane of the sample. 
The sample was placed on a 6-axis positioning stage such 
that during the angular spectroscopy the beam is always 
focused into the same spot. The forward propagating 
THG signal pulses were detected by a photomultiplier 
tube and gate-integrated by an oscilloscope. We used 
the p-p polarization configuration — illuminating with p- 
polarized light and selecting only the p-polarized part 
of forward propagating light before the detector. For 
all measurements spectral filtering (Schott RG610 and 
BG40) before the detector was used for picking up the 
desired wavelength. With these filters the third harmonic 
response was orders of magnitude larger than signals at 
other wavelength, i.e. at the pump wavelength. The aver- 
aged THG signal from the pure Si02 substrate measured 
outside the metamaterial area was at least one order of 



magnitude lower than that from the metamaterial area. 
Contributions from the substrate were therefore safely 
neglected. The principle setup is shown in Fig. [3] 

For numerical simulation an extension of the FMM 
which includes the nonlinear interaction was used. 28 The 
method relies on the undepleted pump approximation 
that ignores the feedback of the nonlinearity induced 
field to the pump field. 29 The approach allows solving 
the problem completely rigorously and permits a reliable 
prediction of the diffracted amplitudes of the third har- 
monic fields. 

The third harmonic intensity was measured and simu- 
lated in the forward zeroth diffraction order with the fun- 
damental wavelength exciting the magnetic resonance. 
The angular spectra of THG are provided in Fig. H(a)- 
(d) for the fundamental wavelengths of 1.49 /im, 1.54 /im, 
1.56 jam and 1.60 /im, respectively. The magnetic res- 
onance position for normal incidence is 1.54 jam. The 
maximum of the THG signal is seen at angles of inci- 
dence around 20°. The appearance of this maximum is 
detailed in the discussion section and is believed to be 
caused by the interference of THG from the individual 
layers forming the fishnet metamaterial. The simulation 
shows an agreement with the experimental values. The 
THG signal is expressed in a pump power-independent 
fashion as derived from the numerical calculations; the 
absolute values of the THG signal are valid only for the 
simulation results while for the experimental data they 
are of the same order of magnitude. The estimation of 
the experimental value of the effective nonlinear suscep- 
tibility is Xmi = 10 _18 m 2 /V 2 which is the same order 
of magnitude than the reference value of bulk gold. 30 



V. DISCUSSION 

Plasmon-enhanced THG at the magnetic resonance of 
fishnet metamaterials was reported previously^ It was 
shown that the THG spectra obey the principles of the 
local- field enhanced nonlinear response. It was proposed 
that the wavelength dispersion of the THG efficiency is 
defined by the spectral line of the magnetic resonance 
cubed. The maximum of THG at the angles of about 
20° can neither be explained by means of dispersion of 
the local field factor at the fundamental frequency, see 
Fig.|4]^e)-(h), nor with the linear transmission character- 
istics at the third harmonic wavelength, see Fig. [3{T)-(1). 
Finally, the position of the maximum does not coincide 
with the angular position of the propagating diffraction 
order appearance as illustrated with the vertical dashed 
lines in Fig. [JJi)-(l). In this section we show that, firstly, 
this feature is caused by retardation effects, and secondly, 
it is specific to the antisymmetric electric current struc- 
ture of the magnetic resonance. 

The observed third harmonic radiation is considered 
to be caused by the nonlinear polarization of gold due 
to anharmonic electron movement. Nonlinearities of 
other substances of the metamaterial are neglected since 
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FIG. 5. Parameters of the model and uncompensated charge 
density distribution in the unit cell of the fishnet metamate- 
rial for (a) antisymmetric and (b) symmetric resonances and 
corresponding far- field radiation patterns. The blue area be- 
tween the gold layers is shown for better understanding of the 
layout, no influence of the dielectric is assumed in the model. 



their x^ -tensor components are several orders of mag- 
nitude smaller than that of bulk gold: Xim(Si02) = 
4.6 • 10- 23 m 2 /V 2 , Xmi(MgO) = 1.0 • 10- 22 m 2 /V 2 and 
Xini(Au) = 7.5-l(T 19 m 2 /V 2 . 30 32 Without further dis- 
cussion of the specific source of that third-order nonlin- 
earity, we describe the motion of electrons of gold at the 
third harmonic wavelength within the conducting layers 
of the metamaterial using a model of weakly coupled os- 
cillators. Within the chosen model the phase difference 
between the oscillators in the two layers dictates whether 
the resonance is antisymmetric — currents in the two lay- 
ers are antiparallel to each other, Fig.[5fa) — or symmet- 
ric — currents are parallel, Fig. [5jb). At the third har- 
monic wavelength this phase difference is assumed to be 
equal to the phase difference of the oscillators at funda- 
mental frequency multiplied by three. For the antisym- 
metric resonance the phase difference is equal to 7r, for 
the symmetric one it is equal to zero (see Appendix A). 

With this knowledge we write down the dynamical 
equations for the charge density at the third harmonic 
frequency: 



3 (r, t) = qo cos3u;£ • 

■\s(z- d -±l)-6 



(1) 



and for the current density: 
j^ s (r, t) = 3uqo sin 3cot- 



(2) 
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d + s 



2 J V 2 
for the antisymmetric resonance and the dynamical equa- 



tions for the charge density: 
p s (r,t) = qocos3uit ■ 
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J d + s\ f d + s 
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(3) 



2 J V 2 
and for the current density: 

jy(T, t) = 3ooqo s'm3ujt • 

~s( z - d -±±) +5 ( z+ d -±^ 



(4) 



for the symmetric resonance. Here 5(y) is the Dirac delta 
function, Q(y) is the Heaviside step function and qo is 
the amplitude of the uncompensated charge oscillations 
at third harmonic frequency. The latter depends on the 
magnitude of the nonlinear polarization and is propor- 
tional to the x^ components and the local field factors 
at the third harmonic frequency Ls ou (6) and fundamen- 
tal frequency L U} (6) cubed. The solution of the potential 
equation 



A-— — 

c 2 dt 2 



A (M) = -Moj(M) 



(5) 



is sought. The problem is considered two-dimensional, 
i.e. ^-independent. Firstly, we consider the antisymmet- 
ric resonance. The solution of Eq. (j5j) could be expressed 
with the retarded potential 



A(r,i) 



Mo 

47T 



dV 



r|/c) 



(6) 



Since H = curlA//io, the magnetic field distribution in 
the far field (r ^> r') is expressed in the cylindrical coor- 
dinates by substitution of Eq. (j2j) into Eq. (j6j) as follows: 



H y = H z 



3uqo sin f3 

IT ' sm 

7rr cos p 

k(d + s) sin f3 
2 

= 0, 



kw y cos f3 



• sin(3ujt — fcr), 



(7) 



where k = 3uo / c and j3 = + tt/2. The angular radiation 
pattern R{&) is defined by the averaged electromagnetic 
intensity which the unit cell of the metamaterial emits 
per unit solid angle as a function of radiation angle. It is 
expressed as follows: 



R(P) 



dP 



r[r ■ [E x H]; 



(8) 



For a plane wave applies r • E x H = vyj iiq/cqH 2 . By 
substitution of Eq. (0 into Eq. (|5J) and time averaging we 
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get the angular radiation pattern for the antisymmetric 
resonance: 



Ras(P) (X 



qo tan f3 sin 



kw y cos j3 



sin 



k(d + s) sin/3 



1 2 



(9) 



The radiation pattern can be evaluated for the symmetric 
resonance in the same way by use of Eq. (j4]) and Eq. (j6j): 



qo tan j3 sin 



kw y cos /? 



• cos 



k(d + s) sin/3 



i 2 



(10) 



The polar plots in Fig. [5] show the normalized angular 
dependences of THG calculated using Eq. ((9]) and Eq. 
([TO]) for the antisymmetric and symmetric resonances, 
respectively, for the same parameters. The dependence 
of qo(0) ~ | (0)| 3 can be expressed for the magnetic 
resonance with a Lorentz spectral line: 



dujo 
00 



(ii) 



relying on the approximation under which the local field 
correction factor is proportional to the absorption con- 
tour function. Then, the central frequency of the reso- 
nance ooo(0) is substituted by the truncated Taylor ex- 
pansion in the form of uo(0) = ujq + Odouo/dO account- 
ing for angular dispersion of the resonance. The angular 
radiation pattern of the third harmonics is straightfor- 
wardly connected to the angular dependence of THG. 
The third harmonic radiation is emitted from each unit 
cell of the metamaterial with the relative phase which 
depends on the angle of incidence of the pump. Radia- 
tion from each cell interferes to compose the diffraction 
pattern. The intensity of each diffraction lobe depends 
on the angle of diffraction via the radiation pattern de- 
pendence. If only the zeroth diffraction order is detected 
then the diffraction angle equals the angle of incidence 
and thus the radiation pattern is probed by measuring 
the angular dependence of THG. Now we use Eq. (|9j) 
and Eq. JTT} to calculate the data on angular-dependent 
THG from the fishnet metamaterial. The function used 
is expressed as follows: 



1(0) = B 



t r / r\\ 1*3 t sv, ■ I Sill 

|L((9)| 3 cot((9)sin ' y 



■ sin 



k(d + s) cos 9 



1 2 



(12) 



The parameters in Eq. [TT] are determined from the lin- 
ear measurements, see Fig. [2j The angular disper- 
sion of the resonance central frequency is duJo/dO ~ 
3- 10 12 rad/(deg-s) and 7 = 0.15 ±0.01 ps _1 (corresponds 



to AA FWHM = 220 nm). The parameter B stands for a cal- 
ibration coefficient that was not measured precisely. For 
w y the SEM- measured value was taken and (d+s) was set 
to 250 nm. The angular dependent third harmonic inten- 
sity function (Eq.[12j) is plotted in Fig.H^a)-(d) with solid 
lines. A good quantitative correspondence is observed 
between the experimental data, the numerically calcu- 
lated data and the modeled dependence. From all the 
parameters only (d + s) differs from the experimentally 
measured one. The main reason is general oversimplifica- 
tion of the model, i.e. not taking the real phase velocity 
of the third harmonic radiation inside the metamaterial 
into account, considering pure symmetric/antisymmetric 
modes, assuming infinitely dense charge and current dis- 
tributions, etc. Nevertheless, the model gives an explicit 
way how one can distinguish between symmetric and an- 
tisymmetric resonances of the metamaterial by means of 
its nonlinear optical response. For the symmetric res- 
onance no local extremum is observed at oblique inci- 
dence whereas the maximum is present in the case of the 
antisymmetric resonance. In terms of effective ten- 
sor components of the metamaterial this means that the 

(3) 

Xyyyy component of the medium at the magnetic reso- 
nance is less pronounced than that at the electric res- 
onance. In correspondence with the general concept of 
metamaterials it makes possible to tailor the relation be- 
tween different tensor components by the proper choice 
of the metamaterial resonance and its parameters. More- 
over, it could be seen from Eq. (9) and Eq. (10) that effec- 
tive nonlinearities of the metamaterial straightforwardly 
depend on its dimensions, namely w y , d, and s in the 
framework of the model. 



VI. CONCLUSIONS 

To conclude, a magnetic resonance contribution to 
third-order optical nonlinearities of the fishnet metama- 
terial was shown. It was achieved by means of measure- 
ments of the third harmonic signal in forward direction 
from a fishnet sample and numerical simulations with a 
nonlinear FMM. Interference of radiation from separated 
third harmonic sources is shown to emerge as a local max- 
imum in the angular spectra of the third harmonic signal 
found at oblique incidence. Antisymmetric oscillations 
of currents, which are the intrinsic property of magnetic 
resonances, are found to be responsible for the particular 
radiation pattern. Based on this an analytical model was 
build. The angular characteristic of the third harmonic 
response from the experiment, the FMM and the ana- 
lytical model were compared. A quantitative correspon- 
dence between these data sets is observed. The results 
contribute to a better understanding of the possibilities 
of the nonlinear properties of optical metamaterials with 
plasmonic resonances of different symmetries. 

The authors acknowledge support from German Re- 
search Foundation (SPP 1391 prority program, NanoGu- 
ide), German Federal Ministry of Education and Re- 
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VII. APPENDIX A 

Here we discuss the phase difference between the 
sources of third harmonic radiation. The source of the ra- 
diation are oscillations in the gold layers at the third har- 
monic frequency. We use a coupled oscillators model with 
a nonlinear extension. Uncompensated charges are in- 
duced at the edges of the thick wires of the metamaterial 
by the external electromagnetic field with a polarization 
along the thin wires as shown in Fig. [5J 33 Charge conser- 
vation implies qi(t) = —qs{t) and q2(t) = —q±(t)- Har- 
monic oscillations of the charge densities in two coupled 
layers can be described as a superposition of two eigen- 
modes of the system — the first one corresponds to codi- 
rectional currents in the layers and the second one cor- 
responds to counterdirectional ones. 34 Consider x\(t) = 
Qi(t) ~ Qs(t) = 2#i(t) for the uncompensated charge at 
the upper fishnet layer and X2(t) = q 2 (t) — q±(t) — 2q2 (t) 
for the lower fishnet layer. The linear dynamics of these 
values is described by the coupled harmonic oscillators 
model: 

x x (t) + 2j Xl (t) + ujIx x (t) + (jx 2 it) = fe iujt (13) 

Here is 7 the damping constant, ujo is the central fre- 
quency of the resonance for an isolated layer, a is the 
coupling constant, / is the oscillator strength and cpo is 
the difference of phases of the exciting fields caused by 
the retardation. The dynamics of the asymmetric mode 
X(t) = x\{t) — X2(t) is described by 

X(t) + 2-rX(t)+u%X(t)-<TX(t) = f(l-e^°)e iut . (15) 

The solution of the equation in the frequency domain is 
expressed as 



X(uj) 



/(l-e^>) 



uJq — uj 2 + 2ijuj — a 



(16) 



In the case when Q-factor of the modes is high enough for 
the condition yja ^> 7 to be held the asymmetric mode 



implies x\{uj) + X2(oj) ~ and arg£i(co>) — arg^2(co') = ir 
as a consequence. 

Now we consider a nonlinear addition to the electron 
movements 

x 1 (t) + 27x 1 (t)+u%x 1 (t) + ax 2 (t) + axl(t) = fe iujt (17) 
and 

x 2 (t) + 2jx 2 (t) + u%x 2 it) + a Xl (t) + ax\ (t) = fe^^ , 

(18) 

where a <C j 2 uJq/ f 2 . This restriction corresponds to the 
experimentally observed low conversion 10 -11 ) from 
the fundamental field to the third harmonic field and 
allows one to use the perturbation theory approach. At 
the magnetic resonance apply x\{t) « — x 2 (t) and only 
one equation have to be considered: 

x 1 (t) + 2jx 1 (t)+ujlx 1 (t)-crx 1 (t) + ax 3 1 (t) = fe iujt . (19) 

The approximate solution is reduced to two terms: 

Xl (t) = x\(u)e lUJt + x[(uj)e l3uJt . (20) 

After substituting the solution into Eq. (fT9|) and calcu- 
lating the multipliers of e lUJt and e l3ut , one gets 



f 



uJq — uj 2 + 2i^uo — a 



(21) 



and 



x[(uj) 



uj 2 — (3(ju) 2 + 6ijuj — a 



Analogously one gets 
x° 2 (oj) = 

and 

x' 2 (u) 



uJq — uj 2 + 2i^uo — a 



(x?H) 3 . (22) 



(23) 



uj 2 — (3cj) 2 + 6iju — a 



(4ny 



(24) 



Since the first multipliers in Eq. (23) and Eq. (25) are 
not resonant and have the same phase, the phase differ- 
ence argo^ — argx 2 is defined by the second multipliers. 
These multipliers are equal to (xJ(o;)) and (^(cj)) for 
the upper and lower layers, respectively. As a conse- 
quence aigXi — argx 2 = 3(argxi(o;) — arg^2(cj)) = 3tt 
which means that at the THG frequency the electrons 
move inside two gold layers out of phase. 
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